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Abstract 

In this paper we review the physical applications of the generalized Lambert func¬ 
tion recently defined by the first author. Among these applications we mention the 
eigenstate anomaly of the Hfj ion, the two dimensional two-body problem in gen¬ 
eral relativity, the stability analysis of delay differential equations and water-wave 
applications. We also point out that the inverse Langevin function is nothing else 
but a specially parametrized generalized Lambert function. 
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1 Introduction 


1.1 The definition of the Lambert W function 


The solutions of the transcendental equation 

xe x = a 

were studied by Euler and by Lambert : see |6]. The inverse of the function on the left hand 
side is called the Lambert function and is denoted by W. Hence the solution is given by W(a). 
If — i < a < 0, there are two real solutions, and thus two real branches of W If we allow 
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complex values of a, we get many solutions, and W has infinitely many complex branches 
pmiHlITE] . These questions are comprehensively discussed by Corless et al. [6]. The question 
why we use the letter “W” is discussed by Hayes HU; see also 0. 

For the remainder of this paper, all parameters will be assumed to be real, and our concern 
will be real-valued functions of real variables. 

The Lambert function appears in many physical and mathematical problems. Some recently 
discussed physical models, however, lead to equations giving exp(x) as a given rational function 
of x : the preceding equation is but a simple special case. In the following section we briefly 
describe some significant applications of W. Then we turn to the generalized Lambert function 
and some applications. 


1.2 Applications of W 


The survey paper of Corless et al. [6] describes a large number of applications of IT. This 
function appears in the combinatorial enumeration of trees, in the jet fuel problem, in enzyme 
kinetics, or in the solution of delay differential equations, just to mention a few areas. Some 
specific applications of W in the investigation of solar wind comes from Cranmer 0, some 
applications in electromagnetic behavior of materials are given by Houari [L3], other appear¬ 
ances in electromagnetics are investigated by Jenn im It is known that Wien’s displacement 
constant can be expressed by IT, and in the discussion of capacitor fields IT also appears [Mlj . 
The Lambert function has applications in quantum statistics, too [33.J. A simple mathematical 
application connects IT to the distribution of primes via the Prime Number Theorem j36j . 


2 The generalized Lambert function 


2.1 Definition 


Now we recall the recent generalization of the Lambert function and then we show its applica¬ 
bility in different physical problems where the classical IT function is insufficient. 

In 2006 T. C. Scott [30] and his co-workers defined the generalized IT function as the solution(s) 
of the equation 

a {x-h)(x - t 2 ) ■ ■ ■ (x - t n ) _ 

(x-s 1 )(x-s 2 )---(x-s m ) 

In [23], where this generalization was studied, this generalized Lambert function is denoted by 

I i r / t 1 t-2 . . . t n \ 

w{ ;a) . 

\Sl S 2 ■■■ S m / 


So, in particular, 

W (;aj — log(a), W (°; a) = W{a), 




2 




w 


t 



t + W (ae - *), 



= s- W 



We now describe various applications requiring the generalized Lambert function. 


2.2 Applications 

2.2.1 A problem in molecular physics 

Corless and his co-workers [6] mentioned that there was an anomaly in molecular physics. 
When physicists tried to calculate the eigenstates of the hydrogen molecular ion (H 2 ), the 
results were not matching with the predictions [10]. The problem was that - being unable to 
solve their equations - the physicists used numerical approximations which were inadequate. 
This problem originally emerged in 1956 and was solved by Scott and his coworkers [27] just in 
1993. In the analytic solution the Lambert function appears which helped to take exponentially 
subdominant terms into account in the solution which could explain the anomaly. Briefly, using 
the double well Dirac delta function model, the wave equation to be solved was the following: 

+ A< ^ X ~~ = E ( X ^' 

The solution then becomes 

^ = Ae~ dlxl + Be~ dlx - Rl . 

Looking for the possible values of d, we make %[} be continuous on each well, and this requirement 
gives that d± must satisfy the transcendental equation 

d± = q{ l±e~ d±R ). 

This can be converted to an equation solvable directly by using W, and we have that 

W(±qRe~ qR ) 


Finally, the eigenenergies of the system are 



This already matches to the predictions. 

If one were to solve a similar eigenvalue problem for ions with unequal charges, W is no longer 
sufficient. Instead, one must solve a transcendental equation of the form 


e cx = a 0 (x - ti)(x - t 2 ). (1) 

The solutions, as it is easy to see, can be given by the generalized Lambert function: 


cW 


ct\ ct 2 c 1 

’ ao, 
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2.2.2 A problem in general relativity 

“One of the oldest and most notoriously vexing problems in gravitational theory is that of 
determining the (self-consistent) motion of N bodies and the resultant metric they collectively 
produce under their mutual gravitational influence... ” [21]. Mann and Ohta were considering 
the equations of motion of two bodies in one spatial dimension (to make the field equations 
tractable and eliminate gravitational waves, this dimensional restriction is very useful). The 
solution contains the W function for the case when the two bodies have equal masses. When 
the masses are different, one needs to solve a more general equation of the form m 

e~ a = ao(z - h)(x - t 2 ). (2) 


Refering back to the previous problem on the eigenstates of the hydrogen ion, we see that we 
have the identical mathematical problem. 


2.2.3 Bose-Fermi mixtures 


Scott mentioned [27] that in some recent physical investigations on Bose-Fermi mixtures the 
transcendental equation 


- c *_ (x - t!)(x - t 2 ) ■ ■ ■ (x - t n ) 

6 — CLq - 


(x - Si)(x -* s 2 ) • • • (x - s m ) 


appears. We see that 

cW | 

gives the solution(s) of equation fl3])- 


ct 1 ct 2 
CSi cs 2 


ct n _ c n ~ m \ 
CSm &0 / 


(3) 


Independently of our study on the generalized W function [23], Scott et al. [281129] studied the 
number of solutions of equations of the form (j3jl . 


2.2.4 Delay differential equations 

The generalized Lambert function helps in the solution of linear constant-coefficient delay 
differential equations (and stability analysis for delay differential equations). 

Consider linear constant-coefficient delay differential equations of the form 

n m 

J2 a kU {k) (t) = J2 b k u {k \t - t). 

k =0 k =0 

Seeking solutions of the form u(t) = exp(A t) leads to the characteristic equation 

n m 

J2a k \ k = e- XT J2 b ^ k - 

k =0 k =0 

Next, suppose the polynomials are factorized to give 

n m 

a n n (^ - 4) = bme~ Xr n(A- S k )• 

k= 1 k =1 
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This is an equation solvable by the generalized Lambert W function. 

To show an application, we cite S. A. Campbell: “Second-order delay differential equations arise 
in a variety of mechanical, or neuro-mechanical systems in which inertia plays an important 
role... Many of these systems are regulated by feedback which depends on the state and/or the 
derivative of the state. In this case the model equations take the form 

u(t) + bu(t ) + au(t) = f(u(t — — r)), 

where a, b are positive constants representing physical attributes of the system, r is the time 
delay, u, u(t — r ) are the values of the regulated variable evaluated at, respectively, times t and 
t — t and the function, f(x,y), describes the feedback.” |3j. Linearizing about an equilibrium 
solution gives a linear constant-coefficient delay differential equations (with n — 2, a 2 = 1 and 
m — 1 in the above) whose characteristic equation is 


(\-t 1 )(\-t 2 ) = b 1 e~ XT (\-s 1 ). 


( 4 ) 


Another application arises in the paper of Chambers [4]. He considered some differential- 
difference equations with characteristic equations solvable by the generalized Lambert function 
with one upper and one lower parameter Efl eq. (2.1b)]. 


2.2.5 Inverting dispersion equations for water waves 

The dispersion relation for periodic water waves propagating into still water is 

to 2 = gk tanh (kh ), 

where to denotes the wave’s frequency, k its wave number, g the acceleration due to gravity 
and h the mean water depth. There are circumstances, e.g. in connection with computations 
using the ‘Mild Slope Equation’, where g, h and uj are given and one needs to solve for k. The 
numerical approximations for doing this are satisfactory, but it seems worth noting that exact 
inversion is possible here, and in more general settings. 

With x = kh and y = uj 2 h/g the preceding dispersion relation is y — xtanh(a;). As a;tanh(x) 
is increasing in x for x > 0, for each y > 0 there is a unique x > 0 solving this equation. 
Some simple bounds on the positive solution x(y) are immediate. As xtanh(x) < min(x, x 2 ), 
x > ma x(y, y/y). The equation rewrites to 

^2x = x + y 
x-y 


( 5 ) 


and its solution is 



(iy /) 

v-% J 


Waves in flows with two layers of different densities have been studied and in 

Pi ~ P2 


p. 421] 


or = gk 


pi coth(fc/i) + p 2 
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Using the same definitions of x and y as in the preceding paragraph the equation above becomes 


y 


= x 


Pi ~ P2 

Pi COtll(x) + p 2 


J 


which rearranges to 


e 2x = 


x + y 


x 


P1+P2 y ' 

pi—P2 y 


This becomes equation (J5J) when p 2 = 0. Once again one can find an expression for x(y) in 
terms of generalized Lambert functions. 


2.2.6 The inverse Langevin function 

The Langevin function L plays an important role in the study of paramagnetic materials and 
of polymers like rubber. More precisely, L describes the magnetization of a paramagnet under 
the presence of outer classical (i.e., non quantical) magnetic field [18] . The expression of L is 


L(x) 


coth(ar)-. 

x 


The inverse of L is widely studied, and good methods of approximating it are available. A 
comprehensive article on the approximations of L~ l is pH] , 


We point out here that the inverse Langevin function L _1 is also a special case of the generalized 
Lambert function. This can easily be seen by the fact that the transcendental equation L(x) = a 
can be rewritten as 


L{x ) 


e 2x + 1 1 
e 2x — 1 x 


After some rearrangement and the variable substitution x —» — \x we get that this equation is 
equivalent to the equation 


e 


— X 


a + 1 X 0+1 
a — 1 x - 

a— 1 


( 6 ) 


Its solution is then given by 


W 


( — 

I a+1 . 

I -2-’ 

\ a—l 


a — 1 
ci T 1 


Hence 


-- 1 , 


a) — —2W | a + 1 ;—! 

a-i 0+1 


or equivalently 

In what follows we collect some information from [2B] on some particular cases of the generalized 
Lambert function. 
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3 The case of one upper and one lower parameter 


Theorem 1 The Taylor series ofW around a = 0 is 


W ( *; a) =t-Tj2 L rc( nT ) e - n V, 


71=1 


n 


where T = t — s 7 ^ 0, and L' n is the derivative of the nth Laguerre polynomial 


Moreover, whent < s the radius of convergence of the power series in Theorem Q] is the following 


r = e* lim 

n—yoo 


4-i m 


L‘i\{n + l)T) 


= ^- 2 ^ 5=3 


3.1 The r-Lambert function 


Mezo and Baricz [23] studied the solution of the equation 

xe x + rx = n. 

This can be rewritten as 


x 


r — e 


x-*’ 

r 


from where it is apparent that the solution is 


W( ° ; —r 
\n/r j 


( 7 ) 


( 8 ) 


This particular (generally multivalued) function is called r-Lambert function by the authors in 
[ 23 ] and was denoted by W r (n). It can be shown that 

lim W I ° ; -r) = W 0 (n ) = Win), 

0 yn/r J 

the classical Lambert function. We also remark that the defining equation (EJ) of W r (n) does 
not contain any singularity around r = 0. 

Depending on the parameter r, the r-Lambert function has one, two or three real branches and 
so the above equations can have one, two or three solutions (we restrict our investigation to 
the real line). 

Theorem 2 Let r 7 ^ — 1 be a fixed real number. Moreover, we define the polynomial Mjf l \y) 


as 


«?>(</) = £ 


k - ( k ' 


n 


-y) 1 - 


2=1 


Then in a neighbourhood of x = 0, 

W r (x) = =r + £ Cl 

r + 1 


^ ^ - r (n) ( 1 


n =2 


X 


r + lj (r + l) n n! 


( 9 ) 


( 10 ) 
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It is known [6] that for the classical Lambert W function 


W(x) rs-/ log(x) + log 


log(a;) 


as x —* oo. The asymptotics of the r-Lambert function is studied in [23j: for any r G 

1 r\ 


W r (x) r^j log(x) + log 


log(x) X 


as x —» oo. Moreover, 


as x —> —oo. 


W r {x) ~ 


1 

-x 

r 


4 The case of two upper parameters 


The Taylor series of W^;oj (see Theorem [T]) contained easily identihcable polynomials: the 
derivative of the Laguerre polynomials. In the case of W ( >2 ; aj we have that the Bessel 

polynomials |19J 

appear (sometimes y n (z) notation is used). The following is a result of Mugnaini 
Theorem 3 The Taylor series ofW t2 \ a^j around a = 0 is 


where T = t 2 — t\. 


Closing remarks 


The first author wrote C code to calculate the real r-Lambert function on all the branches. 
This can be downloaded from https://sites.google.com/site/istvanmezo 81 /. 
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